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Abstract. Recently it has been shown that the zero-energy eigenstate - correspond- 
ing to the stationary state - of a stochastic Hamiltonian with nearest-neighbour 
interaction in the bulk and single-site boundary terms, can always be written in 
the form of a so-called matrix-product state. We generalize this result to stochas- 
tic Hamiltonians with arbitrary, but finite, interaction range. As an application 
two different particle-hopping models with three-site bulk interaction are studied. 
For these models which can be interpreted as cellular automata for traffic flow, we 
present exact solutions for periodic boundary conditions and some suitably chosen 
boundary interactions. 
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1 Introduction 



The so-called Matrix-Product Ansatz (MPA) was developed and used for the con- 
struction of Optimum Ground States (OGS) in quantum spin chains (see e.g. [[I], |2], |3| 
and references therein). Soon the method found its way into the field of one- 
dimensional stochastic processes. The most prominent example in this area is the 
solution of the totally asymmetric simple exclusion process (TASEP) with open 
boundaries by Derrida, Evans, Hakim and Pasquier (DEHP) [|J f], which has reached 
a paradigmatic status in the field of low-dimensional nonequilibrium processes. 
The general principle behind the MPA for stationary states of stochastic Hamilto- 
nians - the so-called cancelling-mechanism - was recognized by Hinrichsen et al ||, 
allowing the treatment of more complex react ion- diffusion systems. The cancelling- 
mechanism may be considered as a generalization of the OGS-concept to stochastic 
systems [0. 

Based on this cancelling-mechanism, Krebs and Sandow (KS) [[| could prove that 
for stochastic processes in one dimension with nearest-neighbour interaction in the 
bulk and boundary fields acting only on the first and the last site, the MPA is not an 
Ansatz, but merely a reformulation of the problem, i.e. the stationary state of such 
processes can always be written as a matrix-product state (MPS). This means that 
- in contrast to the situation for ground states of quantum systems - (generalized) 
OGS are generic for stationary states of stochastic models. 

Many stochastic systems arising in the context of such fields as traffic flow, granu- 
lar matter, chemical reactions and biological motion have naturally an interaction, 
which is not restricted to nearest neighbours. However, the only example of an 
application of the MPA to models with an interaction range r > 2 seems to be the 
work of EBler and Rittenberg ||. These authors introduced a three-site (r = 3) 
version of the cancelling-mechanism used implicitely in M . Preliminary studies per- 
formed by us indicated, however, that this mechanism is not the most general one 
and allows only the solution of some special models. Therefore it was natural to look 
for a generalization of the proposition by KS in order to identify the most general 
mechanism for arbitrary interaction range r > 2. 

Indeed we found that this generalization is possible. In this letter we want to 
show how the mechanism of KS has to be modified to be amenable to systems 
with interaction range r > 2. After clarifying the cancelling-mechanism - which is 
the main result of this letter - the proof of our propositions is a straightforward 
generalization of that of KS. We therefore omit it here and postpone the details to 
a later publication |L(J. Moreover we want to show the usefulness of our method, by 
solving two different particle-hopping models with three-site interaction in the bulk 
for periodic boundary conditions and some suitably chosen boundary interactions. 
Both models are interesting on their own since they can be interpreted as cellular 
automata describing traffic flow on a highway. 



3 At the same time a solution using a different approach has been presented in 0. 
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2 The generalized Krebs-Sandow proposition 



We split the original proposition of Krebs and Sandow in two independent proposi- 
tions A and B. The main reason is that proposition A gives a sufficient local criterion 
for the existence of MPS and covers stochastic processes both with boundary inter- 
actions and periodic boundary conditions. Proposition B deals with boundary inter- 
actions only where it is possible to show explicitly that the criterion of proposition 
A can always be fulfilled. 

Consider a stochastic process on a chain consisting of L sites, where each site can 
be in one of m states. For processes in continuous time the temporal evolution of 
the probability vector is governed by the master equation Q 

d t \P L (t)) = -H L \P L (t)). (1) 
The linear operator is a stochastic Hamiltonian of the form 

L-r+l 

H L (r) = h Mt (r) + J2 hk ,k+l,... ,k+r-l + bright ( r ), (2) 
fc=i 

in the case of boundary interactions. 

For periodic boundary conditions (p.b.c.) it has the form 

L 

= y^fc,fc+l,...,fc+r-l- (3) 
fc=l 

The integer r denotes the range of interaction in the bulk, i.e. hk,k+i,... ,k+r-i ac ts on 
r sites beginning with site k; h\ e f t (r) and /i r i g ht(^) are boundary interactions acting 
on the first, respectively last, r — 1 sites of the chain. Note that the stochastic 
Hamiltonians (Q) and (|3]) are in general non-hermitian. 

From now on we will only be concerned with the steady state solution of (|l|), which 
is the eigenvector \Pl)o of Hl with eigenvalue 0. Since H L is a stochastic matrix, 



at least one such eigenvector exists We will refer to \Pl}o as zero-energy 

eigenvector in the following. 

Moreover we address the question whether \Pl)o can be written as a matrix-product 
state 

\Pl)o = ^-{W\1>° l \V), (4) 
in the presence of boundary interactions. 

For periodic boundary conditions one has to modify the ansatz to guarantee trans- 
lational invar iance: 

\P L ) = ^- Traced]. (5) 



4 For an introduction to master equations in form of an imaginary-time Schrodinger equation, 
see [Ol and references therein. 
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In both cases T> is a vector of dimension m with components Di, where the D{ are 
matrices acting on some auxiliary vector space A. \V) G A and (W\ G A* are vectors 
in A and its dual A*, respectively. 

Proposition A gives a sufficient local criterion for the existence of a matrix-product 
state - the so-called cancelling-mechanism. Before we proceed we have to make the 
following definition: 

Definition. X(r) is a column-vector with m r ~ l entries -X'j 1) j 2) ... ) j r _ 1J , where 2 7 G 
{1,2,... , m}. The X il!i2t ... jir _ 1 are matrices acting on the vector space A defined 
above. The position of X ilii2; _ )ir l in the vector X{r) is given by 1 + {ii — 1) + 
m(i 2 - 1) + . . . + m r ~ 2 (i r _i - 1). 

Using this definition we now can formulate the first part of the generalized KS- 
proposition: 

Proposition (A), (i) If one can find m matrices Di andrrf' 1 matrices X{r) such 
that they fulfil 

r times 
. ^ . 



^fe,fe+i,...,fe+r-iC D ®V®---®V) = X{r) ®V-V® X(r) (6) 

then \Pl)o — Trace[D® L ] is an zero-energy eigenvector of H^\r) (see (^)), i.e. a 
stationary state of the underlying stochastic process. 

(ii) If in addition to ^) one can find vectors \V) G A and (W\ G A* such that 

r—l times 



(W\k(r)(p i V ® • • ■ i V) = -(W\X(r), (7) 
and 

h r {r)( V®V®---®V) \V) =X(r)\V), (8) 



r—l times 



then \Pl)o = (W\V® L \V) is a zero-energy eigenvector of Hl(t) (see ^)) with bound- 
ary interactions hi e f t (r) and bright (^)- 

Proof. The proof is a straightforward generalization of the proof in || . For details 
see 0. □ 



Remark. The relations (Qj and (^) are the most general cancelling-mechanism 
for the stochastic processes considered here. 

In the case of boundary interactions we are able to show more. Here the matrix- 
product state is not an ansatz, but merely a reformulation of the fact that the 
stationary state is a zero-energy eigenvector of Hi, (r ) for all system lengths. 

Proposition (B). Given a stochastic process described by a stochastic Hamiltonian 
of the form ^) which has a unique stationary state for any system length L. Then 
the eigenstate \Pl)o with eigenvalue corresponding to this stationary state can 
be written as a matrix-product state (W\V® L \V) with V = (D 1: D 2 ,... ,D m )* and 
vectors (W\, \V). Moreover one finds m*- 1 * -1 ) matrices X{r), such that the cancelling- 
mechanism is fulfilled. 
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Proof. See [10 . As in [Q one can give an explicit construction for the operators 
involved. □ 



Remark. In the case of m = 2 and r = 3 it is easy to see, that the cancelling- 
mechanism proposed in j^j is a special case of our mechanism. Using this special 
mechanism a proposition similar to our proposition A could be formulated. Nev- 
ertheless our mechanism has the advantage that in case of boundary interactions 
it leads to an algebra which is manifestly non-trivial f\. This stems from the fact 
that in the course of the proof of proposition B we construct an explicit non-trivial 
representation of this algebra. A discussion of this point can be found in JJ^/. 



3 Application of the proposition 

We present now applications of our cancelling-mechanism by determining the sta- 
tionary states of one-dimensional two-state (m = 2) reaction-diffusion systems with 
three-site interactions in the bulk. In the following we investigate two models which 
are generalizations of the TASEP. 

In the first model (model A) particles hop exclusively in one direction, say to the 
right, along a one-dimensional chain of length L with periodic boundary conditions. 
Particles hop one site to the right at rate pi, if this site is not occupied. If a particle 
has two empty sites in front of it, it may also move two sites to the right with rate 
p 2 . The stochastic Hamiltonian of model A has the form ([3]) with r = 3. The local 
operator /ifc,fc+i,fe+2 acts on sites k, k + 1 and k + 2 and its explicit form is given by 
the dynamics of the model. 

Model A is of obvious relevance for the modelling of traffic flow. It can be interpreted 
as a model for cars which have a maximum velocity v max = 2 moving on a single- 
lane highway. Up to now no exact solutions of probabilistic traffic flow models with 
Vmax > 1 are known. 

For the stationary state we make an ansatz of the form (|5|) with T> = (E, D) 1 
where E,D G End(U) and V a vector space. Note that this corresponds to a grand 
canonical description since (|5|) is a superposition of states with different particle 
numbers. However, in our case it is not difficult to obtain results for fixed particle 
densities. 

Using Proposition A from above, \Pl)q is a stationary state if one can find 6 operators 
E, D, Xxx, X\2, X 2 \ and X 22 such that they fulfil the following relations, given by 
the use of our cancelling mechanism 





Pl DEE 
p 2 )DEE 




XuE 


— EXn, 


p 2 DEE 


= X U D 


— EX 12 


X 12 E 


— EX 21 , 


Pl DED 


= X 12 D 


— EX 22 


X 2 \E 


- DX U , 


—piDED 


= X 21 D 


-DX 12 


X 22 E 


- DX 21 , 





= X 22 D 


-DX 22 



(9) 



3 By a trivial algebra we mean an algebra which is equal to 0. 
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The algebra given by 
D,E eR: 



has the following one-dimensional representation with 



D = l-E, 
X 11 = xE/D+p 2 E 2 , 
X 2 i = x - piED, 



X12 = x G K (free parameter), 
X 22 = xD/E- Pl D 2 . 



(10) 



As a consequence the stationary probability distribution is a simple product measure. 
The particle density p is equal to D and the flow J is simply given by 



J{p,Pi,P2) = p(l - p)[pi + 2p 2 (l - p)], 



which gives us the fundamental diagram, i.e. the functional relation between flow 
and density, at hand (see figure [1|). Since model A does not exhibit a particle- hole 

Flow J 




Density p 

Figure 1: Fundamental diagram for model A with p 2 — 1: The dashed line corre- 
sponds to pi = 0.01 and the solid line to p\ = 0.99. 

symmetry, the fundamental diagrams are not symmetric with respect to p = 1/2. 
We have also investigated model A with parallel dynamics, which is more appropriate 
for traffic models [0. Comparing the results of a so-called car-oriented mean- 
field theory (COMF) [14|] with Monte Carlo simulations indicates that COMF gives 
actually the exact fundamental diagram [IU|. 

As a next step we examined model A with open boundaries, i.e. injection of particles 
at the left end of a chain of length L, and removal of particles at the right end. For 
the left end we have chosen the following input rates 



\EE ■ 




DE- 



ED- 



(12) 



\ED ■ ■ ■ 



«3 



\DD ■ ■ ■ 
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The symbol E denotes an empty site and D an occupied one. 
For the right end we have chosen the output rates 



■■■ED\ ^■■■EE\ (13) 




... D D\ P - ^■■■DE\ 

Using the cancelling mechanism one gets the following relations, which have to be 
fulfilled in addition to ([]) 

-(«a + a 2 )EE = -X n , p 1 ED + [3 3 ED = X 11 , 

a 2 EE - a 3 ED = -X 12 , (fc - (3JED = X 12 , 

a.\EE = -X 21 , + fa) ED + (3 A DD = X 21 , 

a 3 ED = -X 22 , -(3 4 DD = X 22 . 



(14) 



We found that the one-dimensional solution ( |10|) serves as a solution of (0)> ^ the 
boundary rates are given by 

ai = v\~D + p 2 DE, p 1 = -p 3 + Pl E + p 2 E(l + E), 

a 2 =p 2 D, f3 2 = -p 3 +Pi +p 2 E, (15) 

a 3 = ati, fa = PiE + p 2 E 2 , 

where f3 3 G R is a free parameter, and 

x =p 1 D 2 E -p 2 DE 3 . (16) 

A more complete investigation of the properties of this model for general values of 
the interaction parameters will be presented in [|TUJ. 

The second model we want to present (model B) has again three-site interactions in 
the bulk. The same model has been studied independently in [TTSJ where equivalent 



results for periodic and open systems have been found. Similar to model A particles 
move along a chain of length L exclusively in one direction. If a particle has two 
empty sites in front of it, it moves one site with rate 1; if only the next site is empty, 
the particle performs the same move with rate A. For A < 1 this model may be 



considered as a traffic flow model with a so-called 'slow-to-start rule' (see |16[ and 
references therein). 

The stochastic process obtained from model B after a particle-hole transformation Q 



is also very interesting flQfl . Here the hopping probability depends on the occupation 



J And a parity transformation so that the particles again move from left to right. 
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number of the site directly behind the particle. If the site to the left is occupied a 
particle moves to the right with rate p%, if it is empty it moves with rate p%. For 
Pi > P2 the particles prefer to stick together which can lead to interesting clustering 
properties. This model might have applications for granular matter and flocking 
behaviour where similar interactions have been studied (see e.g. |p"7f). 
The matrix-product ansatz is the same as for model A and we find for the bulk- 
algebra: 



X n E = EXn, 
DEE = X12E — EX21, 
-DEE = X 21 E - M u , 
X22E = DX21, 



X\\D = EX12, 
XDED = X 12 D- EX22, 
-XDED = X 21 D- DX 12 , 

X22D = DX 2 2- 



;i7) 



For this algebra we found a two-dimensional representation which has a structure 
generic for a 2-cluster approximation [13 



E 



e 1 




D 



d 
1 



ed 



X 



1-X 



X n 
X21 



e 

d 

x 



{d + x)E, 
DE, 



X\2 
X22 



X 



eD + -ED, 
d 

xD. 



(19) 



Remark, (i) It is not hard to prove, that any two operators E and D, which fulfil 



E = eE, 
D 2 = dD, 
XDED = edD, 



(20) 



with e, d being c-numbers and the Xij defined in terms of E, D as in filQ), yield a 
representation of the algebra (17). 

(ii) A first consequence is that any representation of the Temperly-Lieb algebra 
T 3 (±vT) gives a representation of (\P^). This can be seen by defining E = e\ and 
D = pe 2 where e\ and e2 are the generators of T^^y/X) and p is a free parameter. 

With the use of the two-dimensional representation (pl|) it is possible to calculate 
the expectation value of any observable in the stationary state (see e.g. 0). The 
fundamental diagram is given by 



J{p) = P 



1_ v /i_ 4 (l-A)p(l-p) 
2(l-A)(l-p) 



(21) 



As an example we present some fundamental diagrams in figure |2|. Again it is 
non-symmetric due to the lack of a particle-hole symmetry. 
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Flow J 



0.25 
0.2 

0.15 
0.1 

0.05 



* 



♦ ■ 
♦ ■ 



0.2 0.4 0.6 0.8 



Density p 



Figure 2: Fundamental diagram of model B: squares correspond to A 
triangles to A = 0.001 and diamonds to A = 0.5 . 



0.999, 



As for model A, we studied model B with open boundaries, having input rates 

\EE \DE ■ ■ • , (22) 



and output rates 



\ED ■ ■ ■ 



DEI 



■ ■ ■ ED\ 



■■■DDI 



\DD ■ ■ 



— ED\ 



(23) 



■ ■ -EE\ 



■■■DEI 



According to our cancelling-mechanism this gives the following relations, which sup- 
plement the algebra (|X7|) : 



(24) 





(-a x EE\ 








(W\ 


-a 2 ED 
a\EE 
\ a 2 ED ) 




-(W\ 


A12 

X 2 \ 
\X 22 ) 


I 


(3 2 ED 


\ 




}x 11 \ 


fi x DE - f3 2 ED 
foDD - fcDE 


IV) = 


X\ 2 
X 2 \ 




-PsDD 


J 




\x 22 j 



\V). (25) 
Again we found that the two-dimensional representation flTBp, (O) is also a repre- 
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sentation for (pi]), (25) if x is given by 

A 



x± 



1- A 



e + d± \l (e - d) 2 + -ed 



(26) 



where x + and x_ correspond to A > 1 and A < 1, respectively, and the boundary 
rates satisfy 



«i = 1 + 

/?2 = /?3 



.r 



.7' 

ri' 



«2 = A«i, 

— xe 



IV) 



die + x) 

a; 



(27) 



— d — e 



That means that we found a line in the parameter space along which the stationary 
state can be written as a matrix-product state with 2-dimensional matrices. The 
density profiles calculated along this line are flat. 



4 Conclusion 

The generalization of the proposition of KS to stochastic processes with interaction 
range r > 2 leads to the most general cancelling-mechanism for such systems. The 
identification of this mechanism is one of the main results of this letter. Using this 
mechanism the existence of MPS for the stationary state of systems with boundary 
interaction is shown. As an application we were able to solve two models with 
three-site interactions. 

These models are interesting by themselves since they might be of relevance for the 
description of traffic flow or granular matter. We were able to find the stationary 
state of the periodic system which in both cases is given by a finite-dimensional 
representation of the algebra obtained from the MPA. 

For the open system with boundary interactions we could extend the solutions of the 
periodic systems for special values of the input and output rates. For general values 
of the boundary rates one probably needs infinite-dimensional representations. 
Both models have a fundamental diagram with only one maximum. Using the 
argumentation of ll| one can expect that the phase diagrams of model A and B 



with boundary interactions essentially look like the well known phase diagram of 
the TASEP § |. 

A generalization of the KS-proposition to other update schemes J7| has been pre- 
sented in fl9 |. A similar generalization is also possible in our case and will be 



presented elsewhere [ |TU| . 
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